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Abstract 

New flight test maneuvers and analysis techniques for 
evaluating the dynamic response of in-flight thrust models 
during throttle transients have been developed and vali- 
dated. The approach is based on the aircraft and engine 
performance relationship between thrust and drag. Two 
flight test maneuvers, a throttle step and a throttle fre- 
quency sweep, were developed and used in the study. 
Graphical analysis techniques, including a frequency 
domain analysis method, were also developed and evalu- 
ated. They provide quantitative and qualitative results. 
Four thrust calculation methods were used to demonstrate 
and validate the test technique. Right test applications on 
two high-performance aircraft confirmed the test methods 
as valid and accurate. These maneuvers and analysis tech- 
niques were easy to implement and use. Flight test results 
indicate the analysis techniques can identify the combined 
effects of model error and instrumentation response limita- 
tions on the calculated thrust value. The methods devel- 
oped in this report provide an accurate approach for 
evaluating, validating, or comparing thrust calculation 
methods for dynamic flight applications. 

Nomenclature 

A/B afterburner 

AJ nozzle throat area, in 2 

AP area pressure (thrust correlation method) 

2 

Axw aircraft acceleration along its flightpath, ft/sec 

Azw aircraft normal acceleration relative to its 
flightpath, ft/sec 2 

2 2 

A8 nozzle throat area, in or ft 

C D vehicle drag coefficient 
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CIW compressor variable guide vane, deg 
D vehicle total drag, lb 

F excess thrust, lb 

EGT exhaust gas temperature, °C 

EMD engine model derivative 

EPR engine pressure ratio 

g acceleration caused by gravity 

FG gross thrust, lb 

FN net thrust, lb 

FNP net propulsive force, lb 

FN P erroT error in net propulsive force term, lb 

F noz nozzle exterior drag, lb 

FR engine ram drag, lb 

F u inlet spillage drag, lb 

FTIT fan turbine inlet temperature, °F 

FVG fan variable guide vanes, deg 

HIDEC Highly Integrated Digital Electronic Control 
HPVG high-pressure (rotor) variable guide vanes, deg 

IFT in-flight thrust 

L lift, lb 

Mil military power (maximum dry power) 

Max maximum power 

NZ normal load factor, g 

N1 fan rotor speed, rpm 

N2 compressor rotor speed, rpm 

PLA power lever angle, deg 

PLF power required for level flight, deg 

PSO ambient static pressure, IbAn 2 

PS2 engine inlet static pressure, lb/in 2 

PS6 afterburner inlet static pressure, lb/in 2 

PS7 exhaust nozzle inlet static pressure, IbAn 2 

PT4 combustor pressure, IbAn 2 



PT558 

turbine discharge total pressure, lb/in 2 

PT6 

afterburner inlet total pressure, lb/in 2 

RCW 

rear compressor variable vanes, deg 

RTTM 

real-time thrust method 

S 

2 

surface reference area, ft 

SVM 

state variable method (or model) 

777 

inlet total temperature, °R 

TT2 

engine inlet total temperature, °F 

V 

velocity, ft/sec 

w 

aircraft gross weight, lb 

WACC 

digital engine control computed airflow, Ib/sec 

WFT 

fuel flow, lb/hr 

WT 

mass flow temperature (thrust correlation 
method) 

W1 

inlet airflow, lb/sec 

Symbols 

a 

angle of attack, deg 

r 

pitch angle, deg 

(O 

frequency, rad/sec 

p 

density, lb/ft 3 

T 

thrust incidence angle, deg 

Subscripts 

L1H 

left-hand engine (facing forward) 

RIH 

right-hand engine (facing forward) 

t 

time, sec 


Introduction 


This report presents a flight test technique developed for 
evaluating the response of thrust calculation methods dur- 
ing throttle transients. Emphasis is placed on developing 
the calculation methodology and applying the test tech- 
nique. The flight test results presented in this report are 
intended to demonstrate the application and support the 
validity of the test technique. These results are not 
intended to discriminate between any given thrust calcula- 
tion method. 

A primary requirement of aircraft performance testing is 
an accurate measure of in-flight thrust. Various methods 
have been developed to calculate installed thrust indirectly 
from the measurements of related engine parameters. 
Until recently, these methods were designed and therefore 
limited to quasi-steady-state flight conditions at fixed 
throttle settings. New in-flight thrust models have recently 
been developed that promise unsteady flight applications 
and transient throttle capability. Dynamic thrust calcula- 
tion methods require new approaches to in-flight 


evaluation and validation to assess their accuracy during 
dynamic conditions. 

Transient in-flight thrust determination provides new 
opportunities for flight test evaluations. Current aircraft 
performance maneuvers are accomplished purposely slow 
to avoid violating the quasi-steady requirement. 1 Thus, 
today's maneuvers are time consuming and restrictive in 
pilot technique. Improved productivity will result through 
better techniques inherent with the use of transient models, 
albeit an improved understanding of unsteady aerodynam- 
ics will also be required. 2 These developments are essen- 
tial requirements for the vehicle performance evaluation of 
highly maneuverable aircraft. The development of accu- 
rate dynamic simulation models will benefit from the in- 
flight validation of transient thrust models. This develop- 
ment is particularly important in simulators used for such 
pilot-training tasks requiring substantial throttle usage as 
carrier approach and landing, formation flight, and in- 
flight refueling. An accurate transient thrust model also 
provides a means for directly evaluating engine response 
characteristics during throttle or engine control system 
inputs. This ability may be a critical need for developing 
accurate thrust vectoring control systems. 

A number of reports have been written on quasi-steady- 
state, in-flight thrust calculation methods. One of the more 
comprehensive references, In-flight Thrust Determination , 
by the Society of Automotive Engineers (SAE) also 
addresses unsteady influences on thrust determination. 3 In 
general, SAE found that using the traditional thrust calcu- 
lation methods during rapid throttle movement could lead 
to large errors. Such errors are predominately caused by 
transient thermal effects not accounted for in these steady- 
state models. Instrumentation accuracy and response limi- 
tations also influence the dynamic response characteristics 
of the calculated thrust value. Errors in the thrust model 
and in the input variables contribute to the accuracy of the 
thrust value. 

To assess the validity of in-flight thrust calculation dur- 
ing throttle transients, the NASA Dryden Flight Research 
Center, Edwards, California, has developed new flight lest 
analysis techniques based on time-dependent, aircraft and 
engine performance relationships. These techniques were 
applied and evaluated during flight testing of the X-29A 
Advanced Technology Demonstrator (Grumman Aero- 
space Corporation, Bethpage, New York) and the F-15 
Highly Integrated Digital Electronic Controls (HIDEC) 
(McDonnell Douglas Corporation, St. Louis, Missouri) 
aircraft. Although both aircraft were fully equipped for in- 
flight thrust calculation and vehicle performance determi- 
nation, no special attempts were made to optimize instru- 
mentation for the dynamic thrust response. These vehicles 
simply provided a convenient opportunity to apply and 
assess the procedures formulated in this report. 
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This report presents the methodology developed for 
evaluating the dynamic response of in-flight thrust models 
during transient engine operation. Flight test maneuvers 
and analysis techniques are developed and applied to 
high-performance aircraft Right test results for a variety 
of in-flight thrust models, including those used tradition- 
ally for steady-state applications, are used to demonstrate 
and validate the test techniques. 

Methodology 

The methodology developed to evaluate the calculated 
thrust value during a throttle transient is based on the force 
balance relationship between engine net propulsive force, 
FNP , and vehicle drag, D. M Figure 1 shows an aircraft 
force balance diagram. For simplification, the following 
assumptions were made: 

1. All forces pass through the center of gravity. 

2. The aircraft is in symmetrical flight (no sideslip). 

3. The gross thrust, FG, vector has no lateral inclina- 

tion. 

By resolving the forces along the fiightpath (wind coor- 
dinate system, X-axis) and assuming constant mass, the 
following equation results: 

FG • cos (a + x) -FF-F, , 

spil 

-F noz -D-W sin(Y) 

= (Wig) Axw (1) 

Where 

g = acceleration caused by gravity, 
ft/sec 2 

FG - gross thrust, lb 
FR = engine ram drag, lb 
F spil = inlet spillage drag, lb 

F = nozz l e exterior drag, lb 
W = aircraft gross weight, lb 
D = aircraft drag, lb 

Axw = aircraft acceleration along its 
fiightpath, ft/sec 2 

a = angle of attack, deg 
y = pitch angle, deg 
T = thrust incidence angle, deg 

The first two terms in equation (1) also define the net 
thrust, FN , along the fiightpath. That is, 

FN = FG * cos (cx + t) — FR (2) 


Net propulsive force, FNP , is defined as the summation of 
all the propulsion-related forces acting on the vehicle. That 
is, 

FNP = FG • cos (a + t) -FR- F spU - F noz (3) 

The right-hand side of equation (1) defines the aircraft 
excess thrust term, F ex . That is, 

F„ = (Wig) • Axw (4) 

Substituting equations (3) and (4) into equation (1) yields 
F ex = FNP-D-W sin(y) (5) 

Equation (5) provides an important relationship between 
aircraft drag, net propulsive force, gravity component and 
measurable excess thrust of the vehicle along its flight- 
path. The gravity component is relatively small when pitch 
angle is small. Equation (5) also provides a useful method 
for determining vehicle drag from flight data. Because 
excess thrust depends on acceleration measurements, this 
approach for determining vehicle performance is com- 
monly referred to as the “accelerometer method/’ 5 

For the wings level, quasi-steady flight condition, it can 
be assumed that D and W sin(Y) are constant Assuming 
the flight condition remains constant, a perturbation in net 
propulsive force, A FNP, will be directly comparable to a 
change in excess thrust. For example, 

A F ex = A FNP (6) 

Likewise the time-rate-of-change in net propulsive force is 
equal to the time-rate-of-change in excess thrust: 

dF ex fdt = dFNP/dt (7) 

Equations (6) and (7) form the basis for the quantitative 
evaluation of net propulsive force during a throttle tran- 
sient Remember that the net propulsive force parameter is 
comprised of several engine-related forces and not just 
gross thrust. The excess thrust parameter is the reactionary 
response of all the forces acting on the vehicle along its 
fiightpath axis and is directly measurable. 

Equation (6) can be modified to account for a noncon- 
stant drag term if required. For the case where gross thrust 
or net propulsive force does not pass through the center of 
gravity, an additional trim drag term may be required to 
account for the change in elevator trim position required to 
maintain level fiightpath. 
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AF„ = AFNP-AD lnm (8) 

Jet power effects on trim drag can be large when the thrust 
line passes far from the center of gravity, and the change 
in thrust is large. 

If the aircraft changes velocity significantly during the 
maneuver, a time-dependent drag correction similar to 
equation (8) must be used to account for the resulting drag 
change. For example, 

A F ex = &FNP-AD (9) 

Where AD = D ^ ^ - D t = 0 

The time-dependent AD term is determined by first calcu- 
lating the steady-state drag value, D f - o using equation 
(5). Once D, = o is known, the drag coefficient, C D> can be 
computed using the drag equation below. 

D = 1/2 p V 2 S C d (10) 

The time-dependent drag value, D^, can then be deter- 
mined directly from equation (10) using the measured 
velocity and assuming C D remains constant. Two factors 
must be considered when evaluating if the drag coefficient 
is constant: changes in trim angle of attack, a, and com- 
pressibility effects. 

Below approximately Mach 0.7, the primary concern is 
a change in aircraft angle of attack that results in a change 
in the induced drag coefficient. For subsonic operation, the 
induced drag coefficient change can be determined from 
the change in lift coefficient as a function of angle of 
attack. Another factor that must be considered when test- 
ing in the transonic flight regime (above approximately 
Mach 0.7) is compressibility effects or changes in the 
wave drag. The drag coefficient can change significantly if 
the Mach number varies near the transonic drag rise 
region. 

Flight Test Maneuver and Analysis Techniques 

To assess the application of the time-dependent, aircraft 
and engine performance relationships outlined in the 
Methodology section, two flight test analysis techniques 
were developed. One technique is useful in evaluating 
throttle transients or steps. The other technique is useful in 
evaluating the response of a throttle frequency sweep. 
Both techniques are described in the subsections below. 

Throttle Step Maneuver 

The throttle step maneuver and analysis technique is 
based on the engine and aircraft dynamic performance 


relationship presented in equation (6). The maneuver is 
initiated by first stabilizing the aircraft at the desired flight 
condition such that F ex = 0. For single-engine operation, 
the engine is stabilized at the power required for level 
flight, PLF, throttle setting or power lever angle, PLA. 

For aircraft with two or more engines, the test engine 
may be stabilized at any desired initial throttle setting pro- 
vided the remaining engines can be positioned in such a 
way as to achieve an initial steady-state flight condition 
( F = 0). A throttle step input or series of step inputs are 
then performed at various power lever angle rates, ampli- 
tudes, or both, (fig. 2(a)). The change in net thrust during 
the throttle step must equal the change in net propulsive 
force. Large power lever angle inputs will cause the air- 
craft to accelerate and quickly change velocity, thus even- 
tually violating the quasi-steady assumptions (AD = 0). 
One way to minimize this problem is to combine a series 
of opposing throttle step inputs, so an advancing power 
lever angle step is shortly followed by a retarding power 
lever angle step (fig. 2(b)). Even if the quasi-steady 
assumption is violated, equation (9) may still be used to 
evaluate the dynamic thrust value. 

Throttle Step Analysis 

The throttle step analysis method provides a visual com- 
parison of net propulsive force and excess thrust during 
the maneuver. Because of the relationship presented in 
equation (6), a cross-plot of net propulsive force as a func- 
tion of excess thrust will ideally yield a linear relationship 
with a slope of one as illustrated by the dashed line in fig- 
ure 3. Remember that the relationship outlined in equation 
(6) only holds for the short duration where aircraft drag is 
relatively unchanged. Errors in the calculated thrust value 
are depicted by the deviation in net propulsive force from 
the ideal (A FNP = A F ex ) line. When conducting a 
series of power lever angle steps, the curve for excess 
thrust may shift during the fixed throttle portions of the 
maneuver because of changes in vehicle flight condition or 
drag. This shift is illustrated in figure 3. The ideal lines for 
the opposing steps are off-set because of a change in drag 
(and thus excess thrust) because the velocity increases 
while at the maximum power (Max) condition. 

Throttle Frequency Sweep Maneuver 

The throttle frequency sweep method is based on the 
relationships developed in equations (6) and (7). This 
method provides quantitative information on the dynamic 
response of a given thrust model. An approach similar to 
the one used to obtain aircraft stability derivatives in 
flight, where the pilot inputs a range of frequencies to an 
aircraft surface through the stick or rudder, is used. The 
engineer evaluates the aircraft response as a function 
of input frequency. For the throttle frequency sweep 
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maneuver, a range of throttle inputs at increasing or 
decreasing frequencies are used to excite net propulsive 
force which, in turn, affects excess thrust Note that the 
evaluation of thrust model response relative to excess 
thrust does not require any knowledge about the time 
delay between power lever angle and net propulsive force. 

The throttle frequency maneuver is initiated by stabiliz- 
ing the aircraft at a power required for level flight condi- 
tion as outlined in the Throttle Step Maneuver subsection. 
The pilot then conducts a range of constant amplitude 
throttle inputs at various frequencies. Figure 4 shows a 
typical throttle frequency sweep. The input starts at a low 
frequency and progresses to higher frequencies. Con- 
versely, starting at the high frequency and proceeding to 
the low frequency helps to maintain stabilized flight con- 
ditions during the maneuver. 

Throttle Frequency Sweep Analysis 

The throttle frequency sweep analysis requires develop- 
ing a frequency-based relationship between net propulsive 
force and excess thrust. This conversion can be done by 
transforming the data to a frequency domain and con- 
structing a logarithmic plot known as a Bode plot or dia- 
gram. The fast Fourier transformation method is 
commonly used to convert data from a time to a frequency 
domain and was used in this study. The Bode plot (fig. 5) 
provides a direct evaluation of the frequency response 
characteristics (magnitude and phase shift) of a system’s 
output relative to its input as a function of frequency, ox 
For the throttle frequency sweep analysis, the aircraft is 
considered to be the system. When the pilot provides an 
input to power lever angle, the engine responds by chang- 
ing thrust. The aircraft responds with a measurable accel- 
eration. Using excess thrust as the input and net propulsive 
force as the output facilitated obtaining the frequency 
response characteristics of the net propulsive force 
(fig. 5(a)). Excess thrust depends directly on acceleration 
measurements. These measurements are considered to be 
very accurate when compared to the net propulsive force 
calculation. 

The key contributors to errors in the net propulsive 
force value are the model error and input data errors. Input 
data errors can be minimized through proper design of the 
data acquisition system and instrumentation sensors and 
through use of time-correlation procedures to eliminate 
time lags introduced through the data acquisition process. 
The important factor here is that the frequency response of 
the thrust calculation is directly influenced by the fre- 
quency response of its input parameters. Model errors are 
inherent in any analytical representation of a physical sys- 
tem. The Bode plot will show the combined effects of 
model and input errors on the net propulsive force 
calculation. 


The phase value on the Bode plot (fig. 5(b)) quantifies 
the time shift (lag or lead) in degrees between the input 
and output as a function of frequency. For the ideal case 
where the thrust calculation is exactly in phase with the 
measured excess thrust value, the phase shift is 0°. A lag in 
net propulsive force is depicted by a value of less than 0°. 
The magnitude value depicts the amplitude ratio between 
the input and the output as a function of frequency. A deci- 
bel, dB, is equal to 20 Iog 10 of the amplitude ratio. A mag- 
nitude greater than 0 dB represents the case where the 
value of the change in net propulsive force is greater than 
that of the change in excess thrust 

To provide a level of confidence to the frequency 
response results, a data coherence function is calculated to 
evaluate the consistency of the input and output data rela- 
tionship at a given frequency. The data coherence value is 
one when the data exhibits its best correlation of output to 
input This value will reduce when the relationship 
between input and output value is less consistent. If there 
is no dependency between the input and output value, the 
data coherence value is zero. 

One can immediately see the quantitative value of the 
Bode plot for evaluating dynamic thrust calculation meth- 
ods. Such plots provide a method for quantifying the 
dynamic thrust error and can even be used to correct the 
thrust value as a function of the rate of thrust change or 
frequency. Note that these plots provide no insight as to 
the magnitude of the absolute or steady-state thrust calcu- 
lation error and lose their validity at very low frequencies. 

Engine and Aircraft Configurations 

The dynamic thrust analysis techniques were applied 
and evaluated during flight testing of the X-29A Advanced 
Technology Demonstrator and the F- 15 HIDEC aircraft. 
Although both aircraft were fully equipped for in-flight 
thrust calculation and vehicle performance determination, 
no special attempts were made to optimize instrumentation 
for the dynamic thrust response. These vehicles simply 
provided a convenient opportunity to apply and assess the 
procedures formulated here. 

The X-29A aircraft was equipped with an F404-GE-400 
turbofan engine (General Electric Aircraft Engines, Lynn, 
Massachusetts) rated at 16,000 lb thrust, sea level static. 
The engine instrumentation system has previously been 
described. 6 Engine data were recorded at 25 samples/sec. 

The F-15 HIDEC was equipped with an engine model 
derivative (EMD) version of the F100 engine, designated 
the PW1128 (Pratt & Whitney, West Palm Beach, Florida). 
This turbofan engine was rated at 28,000 lb thrust, sea 
level static. A detailed description of the engine and its 
instrumentation system has been given. 7 
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In-Flight Thrust Calculation Techniques 

Four in-flight thrust (IFT) techniques were tested using 
the transient thrust evaluation approaches outlined in this 
report: two traditional gas generator methods, one real- 
time technique, and one dynamic thrust model. All IFT 
models were developed by the respective engine manufac- 
turers or an independent contractor for the specific engine 
and aircraft configurations tested. Their general purpose 
was to provide an accurate calculation of installed gross 
and net thrust for the engine throughout the flight enve- 
lope. Except for the state variable model (SVM), these 
thrust calculation methods were developed and intended 
for steady-state applications only. Their availability and 
widespread use in flight test applications made them con- 
venient methods for this study. 

The thrust values calculated by these models account 
for installation effects as a result of internal performance 
of the inlet, bleed air loss, and shaft power extraction. The 
net propulsive force term accounts for the external forces 
by subtracting the inlet spillage and nozzle drag terms 
from net thrust (eq. (3)). External propulsive drag terms 
are typically less than 1 to 3 percent of the net thrust value. 
Propulsive drag terms depend on aircraft and engine inter- 
actions and vary with power setting and flight conditions. 
Normally these values are determined from wind tunnel 
tests using scaled models with scaled power systems. For 
such cases as the X-29A and F404-GE-400 configuration, 
these values are estimated. 

Manufacturers’ Aerothermodynamic Thrust Models 

The most common approach for determining thrust 
from flight data is to consider the engine as a gas generator 
running at steady-state conditions with thermal and iner- 
tial equilibrium. These so-called gas generator methods 
tend to model the compressor, combustor, and turbine 
components separately to determine mass flow, pressure, 
and temperature conditions at the exhaust nozzle exit 
Internal flowpath measurements within the gas generator 
are used together with mass, momentum, and energy con- 
tinuity principles to calculate flow conditions at various 
stations within the engine and to predict overall engine 
performance. Procedures used to make these calculations 
also use correlations based on steady-state data from 
ground and altitude test facilities. Engine-to-engine varia- 
tions are accounted for by the actual measurement values. 
Successful application of these steady-state correlations to 
flight testing has generally been limited to fixed throttle 
operations at stabilized or quasi-stable flight conditions. 

Two engine manufacturers’ IFT programs based on the 
gas generator method were used in this study. Both the 
General Electric F404 and the Pratt & Whitney F100 
EMD IFT programs use two correlation techniques for 
determining ideal gross thrust: area pressure, AP, and mass 


flow temperature, WT.*’ 9 The thrust value for the area pres- 
sure method strongly depends on an accurate determina- 
tion of nozzle throat area, AS, and nozzle pressure ratio. 
Nozzle pressure ratio is the ratio of nozzle throat total 
pressure to ambient pressure. The mass flow temperature 
procedure requires an accurate determination of engine 
mass flow rate and exhaust gas temperature for its correla- 
tion of thrust and thus requires an accurate afterburner, 
A/B, efficiency model. Table 1 summarizes the input mea- 
surements required for the models presented in this report. 
Detailed development of these gas generator methods for 
gross thrust calculation have been given. 3,10 Gross thrust 
accuracies for these models are expected to be on the order 
of 1 to 2 percent at steady-state conditions for a nominal 
engine. Net propulsive force accuracies are typically 2 to 3 
percent 6 

Real-Time Thrust Method 

A variation of the gas generator method, the real-time 
thrust method (RTTM) was evaluated by NASA Dryden 
on the F404-GE-400 engine-equipped X-29A airplane. 11,12 
Gross thrust is calculated based on a one-dimensional 
isentropic flow analysis in the engine afterburner section 
and exhaust nozzle. The method requires gas-pressure 
measurements from three afterburner locations and a free- 
stream static pressure (table 1). The afterburner pressures 
include the turbine exhaust total pressure and the after- 
burner entrance and exit static pressures. Calibration coef- 
ficients were determined during calibration of the gross 
thrust algorithm from data gathered on the flight test 
engine at the NASA Lewis Research Center (NASA 
Lewis), Cleveland, Ohio, and a ±1.80 percent uninstalled 
gross thrust accuracy was achieved. 

Net thrust is computed from equation (2) where the ram 
drag term is computed from the product of inlet mass flow, 
W1 , and aircraft velocity, V. Inlet mass flow for the RTTM 
was determined by calculating the mass flow rate in the 
afterburner and accounting for bleed losses and fuel flow. 
The RTTM net thrust value was also calibrated against 
NASA Lewis calibration data. An overall ±2.74 percent 
uninstalled net thrust accuracy was achieved. Net propul- 
sive force accuracy is approximately 2 to 4 percent for this 
method. 

State Variable Model 

A SVM combined with aerothermal relations was devel- 
oped for the F100 EMD engine to calculate in-flight thrust 
during steady-state and dynamic engine conditions. This 
model consists of a simplified linear representation of the 
engine rather than the complex representation of each 
component found in nonlinear aerothermodynamic mod- 
els. The SVM is based on a bivariately scheduled piece- 
wise linear, state variable representation of the engine core 
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Table 1. Thrust model measurement inputs. 


X-29A and F404-GE-400 6 



Thrust method 



Area 

Mass flow 

Real-time 

Parameter 

pressure 

temperature 

thrust 

A8 

XX 

X 


EGT 



X 

FVG 

X 

X 


HPVG 

X 

X 


M 

X 

XX 


N1 

X 

X 


PLA 

X 

X 


PSO 

XX 

X 

X 

PS6 



X 

PS7 



XX 

PT558 

XX 

X 

XX 

777 

X 

XX 


WFT 

X 

XX 



F-15 and FI 00 engine model derivatives 


Parameter 


Thrust method 


Area 

pressure 

Mass flow 
temperature 

State 

variable 

AJ 

X 

X 

X 

aw 

X 

X 

X 

M 

X 

X 

X 

N1 

X 

X 

X 

N2 



X 

PLA 

X 

X 

X 

PSO 

X 

X 

X 

PS2 

X 

X 

X 

PT4 



X 

PT6 

X 

X 

X 

RCW 



X 

772 

X 

X 

X 

WACC 



X 

WFT 

X 

X 

X 

TIME 



X 

X - Used in thrust calculation. 


XX - Strongest influence on thrust calculation. 

Note that influence studies have not been performed for the 
F-15 aircraft. 

(fan to turbine). 13 Augmentor and nozzle characteristics 
are modeled with nonlinear aerothermodynamic relations 
which provide input for the conventional mass flow tem- 
perature thrust calculation method. 


The SVM begins by calculating a basepoint data lookup 
and uses a steady-state Newton-Raphson convergence 
algorithm on the first pass to condition input for the 
dynamic portion of the program. Partial derivatives gener- 
ated from relationships to fuel flow and combustor pres- 
sure are used to generate the piecewise linear time- 
dependent thrust calculation variables. 

The manufacturer estimates the steady-state accuracy 
for net thrust to be within 2 to 3 percent and the transient 
accuracy to be within 5 to 10 percent based on model vali- 
dation data. The SVM allows for sample data rates of up to 
50 Hz. Data acquisition was not optimized for this study. 
Table 1 shows the measurement input. Most of these data 
were sampled at 20 Hz. Some engine parameters were 
sampled at 8 Hz because of limitations in the acquisition 
of output parameters from the F100 EMD digital engine 
control system. The low sample rates of the input data will 
obviously limit the response capability of the SVM. 

Throttle Step Results 

Throttle steps at various rates, amplitudes, and direc- 
tions were performed on the F404-GE-400 engine 
installed in the X-29A aircraft beginning at stabilized 
flight conditions. Fast, medium, and slow rates were per- 
formed from power required for level flight to maximum, 
maximum to idle, and military to idle throttle steps. These 
rates varied from about 50 deg/sec for the slow rate to over 
200 deg/sec for the fast rate. Actual throttle rates should be 
adjusted to the specific engine or thrust model being evalu- 
ated. 

Thrust response was evaluated by comparing the change 
in net propulsive force to the change in aircraft-measured 
excess thrust during the transients. Figure 6 shows a sam- 
ple time history of various aircraft and engine parameters 
obtained during a series of throttle transients. Thrust was 
calculated using the mass flow temperature method. 
Remember this method was developed and intended for 
quasi-steady-state applications and is primarily used here 
to illustrate the methodology developed in this report. 

Figure 7 provides a direct comparison of net propulsive 
force and excess thrust for the three throttle rates and 
includes an ideal thrust response line (dashed) based on 
equation (6). Deviations from the ideal response give an 
indication of the errors in the net propulsive force value. 
The results show this thrust method tends to overpredict 
net propulsive force while the throttle is advanced 
(figs. 7(a) and 7(c)) and to underpredict it while throttling 
back (fig. 7(b)). This tendency can be attributed to the 
acceleration scheduling of fuel flow during the advancing 
power lever angle transient. The mass flow temperature 
thrust model uses measured fuel flow, fan inlet tempera- 
ture, and compressor speed to calculate a temperature rise 
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and mass flow rate through the engine. This calculation is 
accomplished using a steady-state fuel-to-air ratio and 
assumes an energy balance between the turbines, fan, and 
compressor components. Because these two assumptions 
do not hold true during a transient, the mass flow tempera- 
ture method overestimates mass flow rate and nozzle 
throat temperature during throttle advancement. Mean- 
while, this method underpredicts both values during throt- 
tle reductions. The tendency becomes increasingly 
pronounced as the throttle transient rate increases. Note 
that the data show closure at the end points (power 
required for level flight, military, maximum, and idle pow- 
ers) once the engine approaches quasi-steady conditions at 
these settings. This finding gives confidence to the validity 
of the throttle step analytical approach. 

Figures 7(a), 7(b), and 7(c) show that a difference 
occurs during the transition between afterburner and dry 
(nonafterburner) power. The primary reason for this differ- 
ence is that the thrust model uses power lever angle to 
select when the afterburner logic is on or off. Because this 
model was intended for steady-state applications, its logic 
does not account for augmentor transition time. This limi- 
tation is most noticeable during the maximum to military 
throttle power step where afterburner fuel flow is prema- 
turely eliminated when the power lever angle reaches mil- 
itary power (fig. 7(b)). As a result, an initial 4000-lb 
reduction in calculated net propulsive force occurs during 
the fast rate transient This reduction amounts to more 
than a 40-percent error from the ideal. 

A similar problem can occur during throttle advance- 
ments from dry to afterburner operation if excessive noise 
exists on the afterburner flow meter when no fuel flows 
though it. In this case, the model uses the noise as the 
afterburner fuel flow value. The result is an overestimation 
of net propulsive force. Figure 7(a) illustrates this problem 
during the power required for level flight to maximum 
power transients. For a short period while afterburner 
operation is being initiated, net propulsive force is influ- 
enced by the incorrect fuel flow value feeding the after- 
burner logic in the model. The net propulsive force value 
takes a sudden drop toward the ideal value as soon as min- 
imum afterburner operation is achieved. This drop occurs 
because the afterburner fuel flow meter is now operating 
normally. 

The dry power step results in figures 7(a) through 7(c) 
are noticeably affected by the engine acceleration and 
deceleration schedule but show little sensitivity to throttle 
rate. The slow rate idle to military power transient 
(fig. 7(c)), differs from the ideal net propulsive force value 
by up to 500 lb. The military to idle power transients 
(fig. 7(b)) deviate by up to 1000 lb. The value for the three 
rate variations is within 100 lb of each other, indicating 


that the throttle may have been at a rate-limiting condition 
imposed by the engine control. Unfortunately because of 
program schedule constraints, a slower throttle rate was 
not performed to confirm this theory. 

Figures 8(a) and 8(b) show a comparison of three thrust 
calculation models obtained using the medium rate throtde 
steps presented in figure 6. In addition to the mass flow 
temperature method, the area pressure method and RTTM 
were evaluated and compared. Bias adjustments were 
made to some methods to improve their agreement at 
steady-state conditions and allow the ideal response line to 
superimpose for each method. 

The area pressure method and RTTM showed improved 
dynamic results over the mass flow temperature method, 
particularly during dry power operation. This result was 
primarily caused by the use of fast-responding pressure 
values for measurement inputs and avoidance of the fuel 
flow measurement The former methods also use power 
level angle as a logical test for afterburner indication. 
Results of such use show the tendency to overpredict 
thrust during throttle advancement from military to maxi- 
mum power and to underpredict thrust while decreasing 
from maximum to military power. Although the area pres- 
sure method and RTTM values were not as accurate as the 
mass flow temperature method during steady-state com- 
parisons, they were better at responding to the dynamic 
step inputs. This finding indicates that a combination of 
methods could yield the improved results for short maneu- 
vers that transition from stabilized to dynamic throttle con- 
ditions. 

Throttle Frequency Sweep Results 

A throttle frequency sweep was performed to evaluate 
throttle rate effects on calculated thrust Data were 
obtained on the F- 15 HIDEC aircraft equipped with two 
F100 EMD engines. The mass flow temperature and area 
pressure gas generator thrust methods used in this evalua- 
tion are similar in theory to those described in the Throttle 
Step Results section. These methods were applied by a dif- 
ferent manufacturer than the results presented in figures 6 
through 8 and use different measured parameters (table 1) 
on an engine with a digital electronic control system. The 
results must be considered independent of those presented 
in the Throttle Step Results section. 

The throttle frequency sweep maneuver began with the 
aircraft stabilized at the specified flight condition. After 
confirming engine stabilization, the left-hand (L/H) engine 
throttle was oscillated at a constant amplitude. The 
oscillation started at a frequency of approximately 0.5 
rad/sec and increased to approximately 8 rad/sec (figs. 9(a) 
and 9(b)). High friction in the throttle lever limited the 
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pilot’s ability to obtain smooth inputs at higher frequen- 
cies. These figures also show how some aircraft and 
engine parameters varied during the maneuver. Figure 9(a) 
reveals that the velocity of the aircraft increased over 50 
ft/sec during the maneuver and resulted in a significant 
change in drag ( D * 0). This velocity deviation was 
compensated for by making an adjustment to aircraft drag 
term as outline in equation (9). 

An equation was developed to calculate the error in net 
propulsive force, FNP error , for the test engine based on 
the ideal relationship presented in equation (5) and 
neglecting the small IV-sin-Cy) term. Note the independent 
treatment of the two engines. 

™ F .rror ( U H, = H + ™ P R,H ~ F „ ~ D <"> 

The aircraft drag and associated drag coefficient were 
calculated during the stabilized portion of the maneuver. A 
time-dependent drag term was then computed from equa- 
tion (10) as a function of the vehicle drag coefficient and 
the velocity squared. The right-hand engine thnotde was 
held constant during the maneuver, and its calculated 
thrust value varied because of the increase in velocity. 
This technique forces the thrust error term to neglect the 
steady-state error of the combined left and right thrust val- 
ues while only determining the dynamic response error of 
net propulsive force. 

Figure 10 shows the resultant error in net propulsive 
force for the mass flow temperature and area pressure 
thrust calculation methods during the frequency sweep 
(fig. 9(a)). The results show a noticeable difference in the 
dynamic response of these two methods. The dynamic 
error in the area pressure method for net propulsive force 
is much more pronounced than that of the mass flow tem- 
perature method. Closer evaluation of the data indicates 
that this error is primarily caused by a time lag in the cal- 
culated net propulsive force area pressure value. The error 
in the mass flow temperature value is much smaller than 
that of the area pressure value. A small phase shift is more 
pronounced at the higher frequencies. 

The approach of using the Bode plot to evaluate the 
thrust response data through the frequency analysis was 
also applied to data from figure 9. Figure 11 illustrates the 
Bode plot results for the area pressure and mass flow tem- 
perature thrust methods. Because of limited stabilized data 
and inherent errors in the Bode plot at low frequencies, the 
data below 0.2 rad/sec should be dismissed. The Bode plot 
results confirm the error analysis results given in figure 10. 
The net propulsive force mass flow temperature value 
stays in phase with excess thrust at frequencies up to 8 
rad/sec. Whereas, the area pressure value indicates a 


significant phase lag (-10°) beginning at a frequency of 
0.5 rad/sec. Confidence in the Bode plot results is greatest 
when the data coherency value is one. A contributor to the 
phase lag in the net propulsive force area pressure results 
may be the engine pressure ratio (EPR) value that is input 
into the model. Pressure measurements are susceptible to 
frequency response errors if these measurements are not 
properly designed for dynamic applications. The EPR 
value shown in figure 9(b) is used as an input to calculate 
net propulsive force area pressure and tends to show some 
time lag when compared to excess thrust Part of the prob- 
lem with the EPR value, PT4/PS0 , is that it is obtained 
from the engines digital electronic control which samples 
this value at a very low rate of 4 samples/sec. Such errors 
provide good examples of how instrumentation response 
error limitations can affect the calculated thrust response. 
Model error can also contribute to such overall results. 

The magnitude results depict the amplitude difference 
between the input and output as a function of frequency. A 
magnitude greater then zero represents when the ampli- 
tude of the net propulsive force value is greater than that of 
the excess thrust value. The Bode plot results show that 
the net propulsive force mass flow temperature magnitude 
values are greater than zero and increase slightly with 
higher frequency. These values indicate a tendency for the 
mass flow temperature method to overshoot during the 
transient. This tendency is probably caused by a model 
error. The engine is on its acceleration and deceleration 
control schedule during the frequency sweep while the 
model assumes steady-state values. Above 8 rad/sec, the 
net propulsive force mass flow temperature magnitude val- 
ues roll off sharply below zero, indicating that the method 
is quickly degrading. 

The area pressure method for obtaining net propulsive 
force results show good magnitude results up to 7 rad/sec. 
At the higher frequency, this method begins to undershoot 
the excess thrust amplitude. In general for this application 
and maneuver, these results show that the mass flow tem- 
perature method is superior in phase shift to the area pres- 
sure method. However, this method shows no 
improvement in magnitude. 

The SVM was evaluated using the throttle frequency 
sweep maneuver. Figure 12 shows the Bode plot results 
from this model. These data were also obtained from the 
maneuver presented in figure 9. Although the SVM results 
show improved magnitude and phase characteristics over 
the area pressure for net propulsive force results (fig. 11), 
they do not show significant improvements over mass flow 
temperature for net propulsive force (fig. 11(b)). The 
coherency value indicates reasonable phase and magnitude 
values were obtained up to approximately 7 rad/sec. 
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Again, limitations in the response of the measured values 
input to the SVM thrust model are a primary suspect in 
these results. Because the flight test analysis method does 
not distinguish between error sources (model or input 
measurements) the expected improved dynamic response 
of the SVM is not reflected in these results. The important 
result here is that the dynamic thrust analysis technique 
can detect limitations in the instrumentation system as 
well as in the model. Once identified, improvements in the 
instrumentation system can be made to improve the 
response of the thrust model. 

As can be seen, the Bode plots provide a comprehensive 
analysis of the dynamic thrust response. The information 
they provide allows for a quantified comparison of thrust 
methods and for the effects of instrumentation system lim- 
itations. Another advantage of the Bode plot is that it can 
be used to correct the dynamic net propulsive force values. 
The quantified magnitude and phase shift errors can be 
applied to the calculated data to provide a more nearly 
accurate value of net propulsive force for an increased fre- 
quency range. The results of this example are limited to 
one flight condition. Additional maneuvers should be used 
to evaluate consistency and changes altitude, Mach num- 
ber, and power level angle amplitude have on these 
results. 

Conclusions 

New flight test maneuvers and analysis techniques for 
evaluating the dynamic response of in-flight thrust models 
during throttle transients have been developed and vali- 
dated. The approach is based on the aircraft and engine 
performance relationship between thrust and drag. Two 
flight test maneuvers, a throttle step and a throttle fre- 
quency sweep, were developed and used in this study, 
Graphical analysis techniques, including a frequency 
domain analysis method, were also developed and evalu- 
ated, Flight test application on two high-performance air- 
craft confirmed the test methods as valid and accurate. 
These maneuvers and analysis techniques were also found 
easy to implement and use. Four thrust calculation meth- 
ods, including those used traditionally for steady-state 
applications, were used to demonstrate and validate the 
test technique. Flight test results indicate the analysis tech- 
niques can identify the combined effects of model error 
and instrumentation response limitations on the calculated 
thrust value. 

Two flight test maneuvers, a throttle step and a throttle 
frequency sweep, were applied to the X-29A Advanced 
Technology Demonstrator and an F-15 Highly Integrated 
Digital Electronic Controls aircraft at the NASA Dryden 
Flight Research Center. Both aircraft were equipped to 


calculate steady-state vehicle performance. These vehicles 
were used to demonstrate the flight test techniques because 
of their availability. The flight results indicate all thrust 
models were affected to some degree by instrumentation 
limitations in frequency response. No attempt was made to 
minimize these errors even though the instrumentation 
was designed for steady-state applications. As expected, 
limitations were found in using the conventional (steady- 
state) thrust models during throttle transients, particularly 
during transitions to and from the afterburner. Model and 
input error limitations were observed. 

Flight test results proved that the throttle step maneuver 
and analysis technique provides a good qualitative evalua- 
tion of dynamic thrust accuracy with limited quantitative 
results. This techniques provides a visual comparison of 
the calculated dynamic thrust value deviation from the 
measured acceleration during a throttle transient. This very 
simple method requires conventional flight test measure- 
ments to implement. Flight test data show good agreement 
between the beginning and ending thrust values and their 
associated excess thrust values and give confidence to the 
validity of the test method. 

The throttle frequency sweep maneuver was demon- 
strated, and the resulting flight data were plotted in the fre- 
quency domain or on a Bode plot. This technique provides 
a complete quantitative evaluation of the dynamic 
response of a thrust model. Although more complex than 
the throttle step method, the results from the Bode plot 
provide more tangible information and a means of correct- 
ing the thrust method phase lag and amplitude variation. 
Results of this effect confirm the flight test techniques as 
valid and provide a useful means for evaluating thrust 
model accuracy during throttle transients. 
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(b) Series of opposing throttle steps. 
Figure 2. The throttle step inputs. 
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Figure 3. Throttle step analysis method. The thrust response error is depicted graphically as a deviation from the ideal. 



Time, sec 


940061 


Figure 4. Typical throttle input during frequency sweep. 
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(b) Bode plot. 

Figure 5. Bode frequency response. Best results for data coherence are obtained when coherence equals one. 
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(a) Aircraft flight conditions. 


Figure 6. Throttle transient time histories for an X-29A airplane with an F404-GE-400 engine. 
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(a) Power required for level flight increasing to maximum power. 
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(b) Maximum to military to idle power. 

Figure 7. Application of the throttle step maneuver to evaluate power level angle rate effects on the dynamic response of a 
thrust calculation model, FNP m for an X-29A aircraft, F404-GE-400 engine, at an altitude of 20,000 ft, and at Mach 0.6. 
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(a) Power required for level flight to maximum power. 

Figure 8. The dynamic response of three thrust models during a throttle step maneuver for an X-29A aircraft, F404-GE 
400 engine, at an altitude of 20,000 ft, and at Mach 0.6. 
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Figure 10. Error in two calculated net propulsive force values determined during a throttle frequency sweep maneuver on 
an F-15 airplane, equipped with an F100 engine ( FNP FNP -D- F a ). 
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Figure 11. Frequency response of net propulsive force to excess thrust for an F-15 airplane equipped with an F100 
(PW1128) engine, at 20,000 ft, and Mach 0.6. 
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Figure 12. Frequency response of the state variable model net propulsive force to excess thrust for an F-15 airplane 
equipped with an F100 (PW 1 128) engine, at an altitude of 20,000 ft, and at Mach 0.6. 
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